arXiv:astro-ph/9808100v2 11 Nov 1999 


Observational Evidence for Superfluidity, and Pinning, in the Core of 

Neutron Stars 


M. Jahan-Miri 

Institute for Advanced Studies in Basic Sciences, Zanjan j5195, IRAN 


ABSTRACT 

The observed large rates of spinning down after glitches in some radio pulsars has 
been previously explained in terms of a long-term spin-up behaviour of a superfluid part 
of the crust of neutron stars. We argue that the suggested mechanism is not viable; 
being inconsistent with the basic requirements for a superfluid spin-up, in addition to its 
quantitative disagreement with the data. One may therefore conclude, for the first time, 
that the presence and the pinned nature of neutron superfluid in the core of neutron 
stars is evidenced by the existing observational data. 
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1. Introduction 

The observational data we are invoking is more 
than a decade old. Also, its implication is rather 
straightforward, once it is realized that there exist 
no alternative explanation for the data, based on the 
role of the (superfluid in the) crust alone. The present 
discussion tends to show the falacy in the only alter¬ 
native possibilty which has been raised to the oppo¬ 
site; thus paving the way to the obvious, though so 
far abondoned, implication of the data on the large 
post-glitch excess spin-down rates of radio pulsars for 
the role of the stellar core. The alternative suggestion 
has invoked a spinning up of the crustal superfluid by 
the spinning down crust (“the container”) over a time 
much larger than the associated relaxation timescale 
(Alpar, Pines & Cheng 1990; hereafter APC90). How¬ 
ever, a closer look at the relative rotation of the super¬ 
fluid and its vortices reveals that the suggested mech¬ 
anism fails quantitatively by, at least, more than an 
order of magnitudes. In addition, the suggested spin- 
up process is also argued to be in contradiction with 
the well-known requirements for a superfluid spin-up. 
This will lead us to the conclusion that the observed 
effect indicates the presence of a superfluid compo¬ 
nent, as well as the pinning of its vortices, in the core 
of neutron stars. We wish to emphasize that the point 
of the present discussion is only to point at the above 
implication of the data, for which no valid observa¬ 
tional evidence has been, so far, invoked (eg., Alpar 
1992). The conclusion of the paper would naturally 
call for, and would be further supported by, a detailed 
and successful modelling of the post-glitch response of 
a pinned superfluid component of the core. The lat¬ 
ter is nevertheless a separate and independent prob¬ 
lem which is not the subject of the present discussion, 
and which has been addressed previously (eg., Mus- 
limov & Tsygan 1985; Chau, Cheng & Ding 1992; 
Sedrakian & Sedrakian 1995; Jahan-Miri 1996). 

2. Overview 

The spin-down rate Q c of the crust of a neutron 
star, with a moment of inertia J c , obeys (Baym et. al. 
1969a) 

I c Cl c = N em - E7if2i (1) 

where N em is the negative electromagentic torque on 
the star, and fi; and A denote the rate of change of 
rotation frequency and the moment of inertia of each 
of the separate components, respectively, which are 


summed over. Steady state implies f2i = fl c = Cl = 
-^p-, for all i, where I = J c +E/i is the total moment of 
inertia of the star. Different models for the post-glitch 
recovery, and in particular the model of vortex-creep, 
aim to explain it in terms of a decoupling-recoupling 
of the superfluid component in the crust of neutron 
stars (Alpar et. al. 1984; Jones 1991; Epstein et. al. 
1992). The rest of the star including the core (su- 
perfluid) is assumed in these models to be rotation- 
ally coupled to the non-superfluid constituents of the 
crust, on timescales much shorter than that resolved 
in a glitch. 

The role of a superfluid component of a neutron 
star in its post-glitch behavior is understood as fol¬ 
lows. Spinning down (up) of a superfluid at a given 
rate is associated with a corresponding rate of out¬ 
ward (inward) radial motion of its vortices. If vortices 
are subject to pinning, as is assumed for the super- 
fluid in the crust of a neutron star, a spin-down (up) 
would require unpinning of the vortices. This may 
be achieved under the influence of a Magnus force 
Fm acting on the vortices, which is given, per unit 
length, as (eg., Sauls 1989) 

F m = -p s K x (v a - i7 L ) (2) 

where p s is the superfluid density, k is the vorticity 
of the vortex line directed along the rotation axis (its 
magnitude k = for the neutron superfluid, where 
m n is the mass of a neutron), and v a and vi, are the lo¬ 
cal superfluid and vortex-line velocities. Thus, if a lag 
u> = — fl c exist between the rotation frequency fl s 

of the superfluid and that of the vortices (pinned and 
co-rotating with the crust) a radially directed Magnus 
force {Fyi)r = PsKvlo would act on the vortices, where 
r is the distance from the rotation axis, and u> > 0 
corresponds to an outward directed (F M ) r , vice-versa. 
The crust-superfluid may therefore follow the steady- 
state spin-down of the star by maintaining a critical 
lag w cr it which will enable the vortices to overcome 
the pinning barriers. The critical lag is thus defined 
through the balancing {Fm)t with the pinning forces. 

At a glitch a sudden increase in Q c would result 
in u> < Ci; cr it, hence the superfluid becomes decoupled 
and could no longer follow the spinning down of the 
star (ie. its container). If, as is suggested (Anderson 
& Itoh 1975), the glitch is due to a sudden outward 
release of some of the pinned vortices the associated 
decrease in Q s would also add to the decreas in to, 
in the same regions. Therefore a fractional increase 
same as the fractional moment of inertia of the 
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decoupled superfluid would be expected. This situa¬ 
tion will however persist only till u> = w cl -it is restored 
again (due to the spinning down of the crust) and 
the superfluid recouples, as illustrated by fi s and f l c 
curves in Fig. 1. The vortex creep model suggests 
(Alpar et. al. 1984) further that a superfluid spin- 
down may be achieved even while lo < w cr it, due to 
the creeping of the vortices via their thermally acti¬ 
vated and/or quantum tunneling motions. A super¬ 
fluid spin-down with a steady-state value of u> < cu cr it, 
and also a post-glitch smooth gradual turn over to the 
complete recoupling for each superfluid layer is thus 
predicted in this model. Notice however that due to 
the varying amplitude of the glitch-induced jump in 
(Alpar et. al. 1984) the sharp recoupling of each 
layer, in the absence of creeping, may occur at a dif¬ 
ferent time. Hence the post-glitch recovery of H c , be¬ 
ing accomplished over an extended time period as the 
different layers of the superfluid recouple, could be¬ 
have smoothly even in the absence of any creep. The 
induced during a superfluid decoupling phase ac- 

i h C 

cording to the creep model is however the same as (or 
slightly smaller than) otherwise. 

3. The Problem 

Recent glitches in Vela, and one in PSR 0355+54, 
have shown values of > 10% with recovery 

timescales ~ 0.4 d, and ~ 44 d, respectively (Lyne 
1987; APC90; Flanagan 1995). The data hence imply 
(Eq. 1) that a part of the star with a fractional mo¬ 
ment of inertia > 10% (- 60%) is decoupled from the 
crust during its observed post-glitch response. This is, 
however, in sharp contradiction with the above glitch 
models, since for the moment of inertia J crust of the 
crustal superfluid Jc y st <: 2.5% (APC90; Chau et. al. 
1993). It should be realized that the disagreement 
with the data is a fundamental shortcoming for the 
crustal models; not just a quantitative mismatch. Be¬ 
cause, the predicted increase for in these models 
is naturally bound to be smaller than the fractional 
moment of inertia of the superfluid (also see the best 
fit results of Alpar et. al. 1993; Chau et. al. 1993); 
except for the possibility raised in APC90 which is 
the point of issue in the following. 

APC90 suggest that the observed large values of 
~ 20% over a time ~ 0.4 d, in Vela, could be ac- 
counted for by assuming that a part of the crustal 
superfluid spins up, over the same period of time. 


The superfluid would then act as a source of an ad¬ 
ditional spin-down torque on the rest of the star and 
could, in principle, result in values much larger 
than its own fractional moment of inertia. For this 
to be realized, a region of the crust-superfluid with 
a moment of inertia / sp and a spin frequency fl sp 
has been assumed to support a tiny (positive) steady- 
state lag w sp = fi sp — fl c ~ 3.5 x 10 -6 rad s , in con¬ 
trast to the much larger steady-state value of the lag 
to > 10~ 2 rad s _1 in the rest of the crust-supcrfluid. 
Hence, a glitch of a size A~ 10 -4 rad s _1 would 
result in a “reversed” situation where > > H sp , and 
which is further suggested (APC90) to be followed by 
a spinning up of the superfluid over a period r sp , as 
indicated in Fig. la. 

4. Quantitative Check 

The above spin-up scenario is however unable to 
account for the observed effect, quantitatively. In 
contrast to the assumption of APC90 that the total 
original frequency difference (fi c — H sp ), induced by 
the glitch, is slowly relaxed during the period r sp we 
argue that only a small fraction of it might be at all 
preserved for any such “long-term” spin-up process. 
This is because the superfluid would be otherwise ro¬ 
tating much slower than its vortices which are, by 
virtue of their assumed pinning, co-rotating with the 
crust (see Fig. la); that is, the superfluid rotational 
lag with its vortices would be much larger than the 
associated critical lag. If so, the pinning could not 
impede the vortex motions (§ 2) and a fast super¬ 
fluid spin-up should be expected. In order to allow 
for the suggested large frequency difference between 
the superfluid and the crust, and at the same time for 
a rotational lag (between the superfluid and its vor¬ 
tices) smaller than the associated critical value one 
needs to further assume that the pinning is “switched 
off”, in contradiction to the assumed pinning condi¬ 
tions. Nevertheless, the superfluid relaxation would 
be still expected to take place very fast, as is further 
discussed below. It may be recalled that the critical 
lag is, by definition, the minimum lag required for 
the Magnus force on vortices to overcome the pinning 
forces. When the instantaneous lag exceeds its criti¬ 
cal value, then the pinning forces (in the azimuthal 
direction) would act further as a major source for 
the torque on the superfluid, resulting in a relaxation 
even faster than in the absence of any pinning (see, 
eg., Tsakadze & Tsakadze 1980; Adams, Cieplak & 
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Glaberson 1985). 

The above argument may be further clarified by 
directly considering the behaviour of the rotation fre¬ 
quency of the vortices (in the region with fi sp ) 
at the glitch, which has not been explicitly specified 
in APC90. We, therefore, discuss the two exclusive 
possibilities that might arise, and could be, in prin¬ 
ciple, physically justified. Both cases lead, however, 
to the same conclusion; intermediate cases for which 
no justification exist should naturally fall in between 
(and there is no indication in APC90 for any spe¬ 
cial effect due to such cases). It may be noted that 
the cases to be considered should not be, however, 
paralleled to the classification of (strong, weak, su- 
perweak) pinning regions, as invoked in the literature 
on the vortex creep model. The latter is based on 
the relative magnitude of the critical lag and reflects 
the long term behaviour of the superfluid relaxation 
toward its steady state lag value. In contrast, the fol¬ 
lowing two cases concerns the instantaneous response 
of the originally pinned vortices upon a sudden jump 
in the rotation rate of the container (see Fig. 1). Ei¬ 
ther, a) the vortices are spun up along with the crust 
(container), that is they remain pinned during the 
sudden spin-up of the container. Or else b) if their 
pinning is temporarily broken, they must relax to a 
state of co-rotating with the superfluid itself. Hence, 
either Hl = H c (pinning conditions), or else Hl = fl sp 
(Hemholtz theorem), just at the beginning of the in¬ 
terval T sp after the jump in fi c (see Fig. la). 

In the case (a) , the superfluid must have been also 
spun up, along with the crust and the vortices, to (at 
least) a state such that H sp — fi c = —w sp (contrast 
with APC90’s picture as depicted in Fig. la); oth¬ 
erwise the pinning would be broken (contrary to the 
assumption) by the associated radial Magnus force 
on the vortices. Experiments on superfluid Helium 
have indeed showed (Tsakadze & Tsakadze 1975; Al- 
par 1978; Tsakadze & Tsakadze 1980) that a pinned 
superfluid either is spun up along with its vessel, or 
it never does so during the subsequent spinning down 
of the vessel (when |fi sp — Hl| < w sp ). Also accord¬ 
ing to the vortex creep model (see Eq. 6 in Cheng 
et. al. 1988), the superfluid spin-up (-down) timescale 
should be extremely short (^ 10”' s for the above 
case in Vela) given that |fi sp — Hl| > w sp , as would 
be the case if only the vortices, but not the superfluid, 
are assumed to be spun up with the crust. 

The case (b), on the other hand, is not in accord 
with the general pinning conditions assumed in the 


vortex creep model, and is not likely to be invoked 
in that context. Nevertheless, the superfluid spin- 
up in the crust of a neutron star, for such a case of 
free unpinned vortices, is again expected to occur over 
very short timescales. The longest timescale for the 
spin-up of the crust by freely moving vortices, due 
to nuclear scattering alone, has been estimated (Ep¬ 
stein et. al. 1992) to be only < 5 s, for the Vela 
pulsar. Such a fast spin-up (-down) of the superfluid 
in the crust of a neutron star, by the freely moving 
vortices, is in fact invoked in the vortex creep model 
as the cause of glitch itself (Anderson & Itoh 1975). 
Therefore, and in either cases (a or b) , the superfluid 
would be spun up to (at least) a frequency such that 
fi sp — fi c = — Wsp, within a period of only few sec¬ 
onds while the jump in H c takes place at the glitch 
(see Fig. lb). Notice that the steady-state lag w sp ac¬ 
cording to the vortex creep model would be slightly 
smaller than the critical lag (though slightly larger in 
the absence of any creeping). The difference is how¬ 
ever a tiny fraction of the critical value itself (Alpar 
et. al. 1984), and is therefore neglected in the present 
discussion. 

Hence, the frequency difference AH between the 
crust and the superfluid, at the start of r sp , would be 
AH = w sp ~ 3.5 x 10” 6 rad s” 1 . This is the total 
available range of rotation frequency that one might, 
in principle, assume to be equilibrated any further 
between the crust and the superfluid. In contrast, a 
value of AH = AH C ~ 1.3 x 10” 4 rad s” 1 has been 
used in APC90. Assuming, for the time, that any fur¬ 
ther spin-up of the superfluid may occur, its duration, 
r sp , may be estimated from (Bayrn et. al. 1969a; see 
also Eq. 2b in APC90) 



I sP AH 

I 'Tsp 


( 3 ) 


where (AH c ) sp is the magnitude of the change in H c 
due to the spinning up of the superfluid. Adopt¬ 
ing the same parameter values as in APC90, ie. 
(AH c ) sp /H c = 0.2, ^ = 5.3 x 10” 3 , AH = w sp = 
3.5 x 10” 6 rad s” 1 , and H c = 9.5 x 10” 11 rad s” 2 for 
Vela, one derives from Eq. 3 a value of 


r sp < 0.3 hr, (4) 

which is too short as compared with the observed 
timescales ~ 0.4 d. Thus, the crust-superfluid can¬ 
not be the cause for the observed large spin-down 
rates even in the case of 1988 glitch of the Vela pul¬ 
sar, addressed in APC90. The disagreement between 
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the predicted and observed timesacles would be even 
worse for the case of 1991 glitch of the same pul¬ 
sar, having observed values of A Cl c /tl c ~ 60% over 
a similar relaxation time (Flanagan 1995). Also, an 
attempt to apply the same crust-superfluid spin-up 
scenario to the case of PSR 0355+54 would result in 
more than three orders of magnitudes difference be¬ 
tween the predicted r sp and the observed relaxation 
time ~ 40 days. 

5. Superfluid Spin-up 

Moreover, the suggested spin-up scenario of APC90 
should be dismissed at once since the required torque 
on the superfluid, during r sp , may not be realized at 
all, under the assumed conditions (see Fig. la). That 
is the pinned superfluid could not be spun up by the 
crust (ie. its container) while the latter is spinning 
down. This is simply because a spinning down ves¬ 
sel (or even one with a stationary constant rotation 
rate) albeit rotating faster than its contained super¬ 
fluid could not result in any further spin up of the vor¬ 
tex lattice which is, by virtue of the assumed pinning, 
already co-rotating with it!. As is well-known, an in¬ 
ward radial motion of the vortices, associated with a 
spin-up of the superfluid, may only occur in the pres¬ 
ence of a corresponding forward azimuthal external 
force acting on them. This is indeed a trivial fact, con¬ 
sidering that any torque on the bulk superfluid may 
be applied only by means of the vortices. However, 
no forward azimuthal force (via scattering processes 
between the constituents particles of the vortex-cores 
and the crust) may be exerted by the spinning down 
crust on the vortex lattice which is already co-rotating 
with it. Thus, no superfluid spin-up might be ex¬ 
pected to occur during the interval r sp , namely after 
the spinning up of the crust and vortices has been 
accomplished (compare Fig. la with Fig. lb). Notice 
that the requirement for a superfluid spin-up, namely 
that the vortices be rotating slower than the crust 
while the latter is itself slowing down, is inherently 
inconsistent with the assumed pinning conditions. 

The suggested long-term spin-up in APC90 is a 
generalization of the vortex creep model for the case 
of a negative lag, in contrast to the usual post-glitch 
spinning-down due to a positive lag invoked in the 
model (eg., Pines & Alpar 1992). However, according 
to the existing formulation of the vortex creep model, 
a spinning up of the superfluid would require a pos¬ 
itive torque N em acting on the whole star (see, eg., 


Eqs. 28, and 38 in Alpar et. al 1984). Application of 
the same formalism to the suggested spin-up episode 
in presence of the given negative N em (which is at¬ 
tempted in Eq. 5 of APC90) is not, a priori, justified; 
it is contradictory. 

The vortex creep model suggests that a radial Mag¬ 
nus force, due to a superfluid rotational lag, results in 
a radial bias in the otherwise randomly directed creep¬ 
ing of the vortices (Alpar et. al. 1984). This might 
be, mistakenly, interpreted to require that given a 
negative lag the inward creeping motion of the vor- 
tice, hence a superfluid spin-up, should necessarily 
follow, irrespective of the presence or absence of the 
needed torque on the superfluid, ft should be how¬ 
ever noted that the role of driving the vortices in¬ 
ward, ie. spinning up of the superfluid, may not be 
assigned to the Magnus force. The Magnus force as¬ 
sociated with the rotational lag is a radial force and 
is also exerted by the superfluid itself, both proper¬ 
ties disqualifying it for being the source of a torque. 
Thus, the point which is to be emphasized here is 
that the obvious requirement for a spin-up process, 
namely the corresponding torque, is indeed missing in 
the suggested mechanism of APC90. Accordingly, the 
inability of the superfluid to be spun up, in this case, 
should not be understood as a direct consequence of 
the pinning preventing a vortex radial motion, which 
could be possibly cured by any assumed radial creep¬ 
ing motions. Rather, the vortices under the assumed 
conditions do not have any “tendency” for an inward 
radial motion, in spite of the presence of an inward 
radial Magnus force which is balanced by the pinning 
forces. A change in the spin frequency of a superfluid 
involves not only a radial motion of the vortices but 
also a corresponding azimuthal one, as may be also 
verified from the solution of equation of motion of 
vortices during a superfluid rotational relaxation (see 
Eq. 9 in Alpar & Sauls 1988 and Eq. 4 in Jahan-Miri 
1998). The torque may be transmitted during only 
such an azimuthal motion and would nevertheless re¬ 
quire and initiate a radial motion as well. Therefore, 
purely radial creeping of the vortices, which is im¬ 
plied by the existing formulation of the vortex creep 
model, may not be invoked as a spinning-up mech¬ 
anism during the transition from H sp — Hl = —w sp 
to H sp ;> Hl- The superfluid would rather remain 
decoupled at a constant value of H sp (if not spinning 
down) during this transition which is achieved due 
only to the spinning down of the crust, as illustrated 
in Fig. lb. 
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6. Concluding 

Decoupling of (a part of) the moment of inertia 
of the crust at a glitch, from the rest of the corotat¬ 
ing star, could readily account for the excess post¬ 
glitch spin-down rates comparable to the fractional 
moment of inertia of the decoupled part. The same 
preliminary fact applies to a (partial) decoupling of a 
(pinned) superfluid component in the crust, or else¬ 
where in the star, for that matter, as well. As sug¬ 
gested by APC90, a component (say, in the crust) 
may result in an even larger excess spin-down rate 
for the star if the decoupled component is further as¬ 
sumed to be spinning up while the rest of the star is 
spinning down; ie. a negative coupling instead of a 
mere decoupling. Nevertheless, as we have shown the 
suggested mechanism of APC90 for such a negative 
coupling of a pinned superfluid part in the crust of a 
neutron star not only fails quantitatively to account 
for the observed effect in pulsars, it is also ruled out 
conceptually since the required torque on the super¬ 
fluid could not be realized at all. 

Thus, the (pinned) crust-superfluid may not be in¬ 
voked as the cause of the observed large changes in 
the spin-down rates of neutron stars. One is therefore 
left to conclude definitely that the observed effect, 
over timescales of a day and more, is caused by the 
stellar core, having a large moment of inertia / core , 
with ~ 90 %. This is, however, by itself an ob¬ 
vious observational evidence for the supej'fluidity and 
at the same time pinning of the vortices in the core 
of a neutron star, given the current “standard” pic¬ 
ture of the interior of these stars (Sauls 1989; Pines 
& Alpar 1992). For, a non-superfluid core would cou¬ 
ple to the crust on very short timescales (< 10 -11 s) 
(Baym et. al 1996b), and could not have any footprint 
left in the observed post-glitch relaxation. Also, a 
core-superfluid with free (unpinned) vortices would be 
again expected to have very short coupling timescales 
of the order of less than one or two minutes (Alpar & 
Sauls 1988; Pines & Alpar 1992; Jahan-Miri 1998). 

The above conclusion, for the presence of pinned 
core superfluid, being a requirement of the observa¬ 
tional data is obviously independent of any specific 
model for the pinning in the core and the post-glitch 
effects of the superfluid relaxation. It is appreciated 
that if the core-superfluid is assumed to be subject to 
a pinning as a whole and also to be completely decou¬ 
pled at a glitch then the induced post-glitch values of 
would be much larger than what has been usu- 
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ally observed for the glitching pulsars. Nevertheless, 
there is no requirement for the pinning of the core 
superfluid that it has to comply to such tight restric¬ 
tions. Given the large moment of inertia of the core, a 
pinning of some parts and/or a partial decoupling of 
the pinned component may as well be invoked so that 
it would result in post-glitch spin-down rates similar 
to that observed. Indeed, a pinning between the vor¬ 
tices and the proton fluxoids (Muslimov & Tsygan 
1985; Sauls 1989) might offer a resolution, as is in¬ 
dicated by our preliminary results (Jahan-Miri 1996; 
Jahan-Miri 1999), even if all of the vortices in the 
core are assumed to be involved. The novel feature 
of this pinning mechanism, as it affects the rotational 
dynamics of the core superfluid, is the moving nature 
of the pinning sites which might allow for a spinning 
down of the superfluid even while the vortices are kept 
pinned! Also, the model of “proton-vortex cluster” 
(Sedrakian & Sedrakian 1995), invoking a somewhat 
different picture of the pinning in the core of neutron 
stars, might serve as another candidate mechanism. 
In any case, while the effects due to the pinned core- 
superfluid on the observable post-glitch relaxation of 
radio pulsars needs to be further demonstrated, the 
observed large excess spin-down rates is pointing at 
an apparent role played by the stellar core, for which 
it ought to be in a pinned superfluid state. 

Raman Research Institute is acknowledged for kindly 
providing financial support as well as the research en¬ 
vironment where this work could be initiated. 
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Fig. 1.— As is shown, schematically, by the thin and 
the dotted lines, a superfluid, as in the crust of a neu¬ 
tron star, whose vortices are pinned and co-rotating, 
at a rate O c , with its spinning down vessel, may spin- 
down at the same rate as the vessel provided a cer¬ 
tain lag between its rotation frequency fl s and fl c is 
maintained. A decrease in the lag, say at a glitch, 
results in a decoupling of the superfluid and hence a 
larger spin-down rate for the vessel. The superfluid 
starts spinning down again once the lag is increased to 
its steady-state value. The thick line represents the 
superimposed behavior of a second superfluid com¬ 
ponent, fl sp , which supports a much smaller steady- 
state lag. Its behaviour as implied in APC90 is shown 
in the panel (a), in contrast to that suggested here, 
shown in (b). The unresolved rising time of fl c at the 
glitch instant t = t g should be reckoned as a measure 
of the expected coupling time of the superfluid due to 
free motion of vortices, in contrast to the suggested 
spin-up period r sp in (a). 
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